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ABSTRACT 

We report the detection of CO J=2— >1, 5— s>4, and 6— >-5 emission in the highest-redshift submil- 
limeter galaxy (SMG) AzTEC-3 at z=5.298, using the Expanded Very Large Array and the Plateau 
de Bure Interferometer. These observations ultimately confirm the redshift, making AzTEC-3 the 
most submillimeter-luminous galaxy in a massive z~5.3 protocluster structure in the COSMOS 
field. The strength of the CO line emission reveals a large molecular gas reservoir with a mass of 
5.3x 10 10 (aco/0.8) Mq, which can maintain the intense 1800MQyr _1 starburst in this system for 
at least 30 Myr, increasing the stellar mass by up to a factor of six in the process. This gas mass 
is comparable to 'typical' z^2 SMGs, and constitutes >80% of the baryonic mass (gas+stars) and 
30%-80% of the total (dynamical) mass in this galaxy. The molecular gas reservoir has a radius of 
<4kpc and likely consists of a 'diffuse', low-excitation component, containing (at least) 1/3 of the 
gas mass (depending on the relative conversion factor aco) ; and a 'dense', high-excitation compo- 
nent, containing ^2/3 of the mass. The likely presence of a substantial diffuse component besides 
highly-excited gas suggests different properties between the star-forming environments in z>4 SMGs 
and z>4 quasar host galaxies, which perhaps trace different evolutionary stages. The discovery of a 
massive, metal-enriched gas reservoir in a SMG at the heart of a large z—5.3 protocluster considerably 
enhances our understanding of early massive galaxy formation, pushing back to a cosmic epoch where 
the Universe was less than 1/12 of its present age. 

Subject headings: galaxies: active — galaxies: starburst — galaxies: formation — galaxies: high- 
redshift — cosmology: observations — radio lines: galaxies 



1. INTRODUCTION 

Our understanding of the physical properties of sub- 
millim eter galaxies (SMGs; see review by Blain et al. 
120021 ) is of key importance to studies of the early for- 
mation and evolution of massive galaxies, as they are 
the likely progenitors of the most massive galaxies in the 
present-day universe. SMGs typically represent compact, 
intense (>1000M Q yr~ 4 ), rather short-lived (<100Myr) 
starbursts with rapid gas consumption through high 
star formation efficiencies that are commonly associated 
with ongoing major mergers. Their star formation rates 
(SFRs) exceed those of 'normal' high-z galaxies with 
comparable stellar mass (-M*) by m ore th an an order of 
magnitude at z^2 (e.g., Daddi et al. 120091 ) . making them 
a comparatively rare, but cosmologically important pop- 
ulation. The extreme star formation events in SMGs are 
typically associated with large amounts of gas and dust 
which often obscure the stellar light and star formation at 
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rest-frame ultraviolet (UV) to optical wavelengths, mak- 
ing their identification at such wavelengths notoriously 
difficult. 

The most insightful way to study SMGs and their star 
formation properties thus usually is through the dust- 
reprocessed UV light from newly formed stars that is re- 
emitted in the far-infrared (FIR) continuum (a measure 
of the SFR), and through line emission from molecular 
gas (typically CO), the fuel for star fomation. Molecular 
gas was detected in >20 SMGs to date, revealing large 
gas reservoirs o f >10 10 Mq in most cases (see Solomon 
& Vanden Bout 120051 for a r eview ). 

Recently, Capak etal. (|2010j ) discovered AzTEC- 
3, an SMG at an unprecedented redshift of z—5.3. 
AzTEC-3 is not only the most distant SMG known 
to date, but also resides in a massive, overdense, 
proto-cluster environment extending out to >2Mpc, 
with two companions within 12.2 kpc. AzTEC-3 
hosts a 1800 M yr -1 starburslQ exhibiting a FIR lu- 
minosity of Lfir=(1.7±0.8)x10 13 L q that likely pro- 
duced a substantial fraction of its current stellar mass 
of M*=(l.O±O.2)xlO lo M . The unusual nature of 
AzTEC-3 and its potential major implications for mas- 
sive galaxy formation at very early cosmic epochs has 
initiated an in-depth study of its physical properties and 
cosmic environment. 

In this Letter, we report the detection of CO( J=2— >1), 
CO( J=5-s-4) and CO( J=6->5) emission toward AzTEC- 
3 (z=5.298), using the Expanded Very Large Array 
(EVLA) and the Plateau de Bure Interferometer (PdBI) . 
We use a concordance, flat ACDM cosmology through- 

8 Assuming a Chabrier (2003) stellar initial mass function (IMF). 
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Fig. 1.— Velocity-integrated EVLA/PdBI maps of the CO(J=2->l) (left), CO(J=5^4) (middle) and CO(J=6^5) (right) line emission 
over 950/852/874 kms -1 toward AzTEC-3. At resolutions of 3.0"x2.2", 3.8"x2.7" and 3.3"x2.1" (as indicated in the bottom left corner 
of each panel), the emission remains unresolved. The cross indicates the position of the sub-millimeter continuum emission (Younger etal. 
120071 ). Contours are shown in steps of lcr=0. 045/0. 10/0. 22 mjy beam" 1 , starting at ±2cr. 



out, with #0= 71 kms 1 M pc 
(Spergel et al. I200H 120071) . 
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2. OBSERVATIONS 
2.1. Plateau de Bure Interferometer 

We observed the CO(J=5^4) (z/ rest = 
576.2679305GHz), and CO(J=6^5) (691.4730763GHz) 
emission lines toward AzTEC-3 using the PdBI. At 
z=5.298, these lines are redshifted to 91.5001 and 
109.7925GHz (3.3 and 2.7mm). Observations were car- 
ried out under good 3 mm weather conditions in the 6C 
configuration on 2010 April 1 and 2, resulting in 4.4 and 
3.1 hr on-source time for the CO J=5— s>4 and 6— >5 lines, 
respectively. The nearby quasar B0906+015 (distance 
to AzTEC-3: 12.8°) was observed every 22.5 minutes for 
pointing, secondary amplitude and phase calibration. 
For primary flux calibration, the standard calibrators 
MWC349 and 3C84 were observed, leading to a cali- 
bration that is accurate within <10%. Observations 
were set up using a total bandwidth of 1 GHz (dual 
polarization; corresponding to —3300/2700 kms -1 at 
3.3/2.7 mm) with the current correlator, and a total 
bandwidth of 3.6 GHz (dual polarization) with the new 
WideX correlator (recorded simultaneously). 

For data reduction and analysis, the GILDAS pack- 
age was used. All data were mapped using 'natural' 
weighting. The CO J=5— >4/6— !>5 data result in a fi- 
nal rms of 0. 52/1. 27mJy beam -1 per 33/27 kms -1 
(10MHz) channel, and 0.05/0. 13mJy beam" 1 
(0.028/0.066mJybeam -1 ) over the entire 1 GHz 
(3.6 GHz) bandwidth. Maps of the velocity-integrated 
CO J—5~ >4/6->5 line emission yield synthesized clean 
beam sizes of 3.8"x2.7" and 3.3"x2.1" and rms noise 
values of 0.10/0.22mJybeam -1 over 852/874 kms -1 
(260/320 MHz). 

2.2. Expanded Very Large Array 

We observed the CO(J=2^1) (iw = 230.53799 GHz) 
emission line toward AzTEC-3 using the EVLA. At 
z=5.298, this line is redshifted to 36.6049 GHz (8.2 mm). 
Observations were carried out under good 9 mm weather 
conditions in D array on 2010 May 24 and 30, result- 



ing in 6.3 hr on-source time with 16 antennas (equivalent 
to 2.2 hr with 27 antennas) after rejection of bad data. 
The nearby quasar J1018+0530 (distance to AzTEC- 
3: 5.4°) was observed every 7 minutes for pointing, sec- 
ondary amplitude and phase calibration. For primary 
flux calibration, the standard calibrator 3C286 was ob- 
served, leading to a calibration that is accurate within 
<10%. Observations were set up using a total band- 
width of 256 MHz (dual polarization; corresponding to 
-2100 kms" 1 at 8.2mm) with the WIDAR correlator. 

For data reduction and analysis, the AIPS package was 
used. All data were mapped using 'natural' weighting. 
The data result in a final rms of 0.20 mjy beam -1 per 
49 kms -1 (6MHz) channel, and 0.030 mjy beam -1 over 
the entire 256 MHz bandwidth. Maps of the velocity- 
integrated CO J=2— >1 line emission yield a synthesized 
clean beam size of 3.0" x 2.2" at an rms noise level of 
0.045 mjy beam -1 over 950 kms -1 (116 MHz). 

3. RESULTS 

We have detected CO(J=2-H), CO(J=5-^4) and 
CO(J=6->5) line emission toward the z=5.298 SMG 
AzTEC-3 at 8a, 10a, and la significance (Fig. [JJ). The 
CO J=5— >4 and 6— >h data yield a combined significance 
of 12<7. We do not detect the underlying continuum emis- 
sion down to la limits of 0.13/0.31 mjy beam -1 (WideX 
correlator: 0. 08/0. 14mJy beam -1 ) at 3. 3/2. 7mm (rest- 
frame 520/434 /im). A 2-3cr peak is present in the 
8.2 mm map (rest-frame 1.3 mm), but offset from the CO 
peak position by ~1". We thus consider the underly- 
ing continuum at this wavelength undetected down to 
-0.09 mjy beam -1 . 

From Gaussian fitting to the line profiles (Fig. [5]) , we 
obtain CO J=2— !>1, 5— >4 and 6— >5 line peak strengths of 
O.45±O.O70 1.78±0.17 and 2.64±0.36mJy at a FWHM 
of 487±58 kms -1 , centered at a (weighted median) 
redshift of z=5.2979±0.00 04 (c onsistent with the opti- 
cal redshift; Capak etal. 2010). This corresponds to 
velocity-integrated emission line strengths of 0.23±0.03, 
0.92±0.09, and 1.36±0.19Jy kms -1 , i.e., line luminosi- 

9 The fit also suggests an underlying continuum component of 
0.09±0.04 mjy, consistent with the marginal peak seen in the map. 
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Fig. 2.— Left: EVLA/PdBI CO(J=2-s-l) (top), CO(J=5->4) {middle) and CO(J=6->-5) (bottom) spectra of AzTEC-3 at 6/10/10 MHz 
(49/33/27 kms" 1 ) resolution (histograms), along with Gaussian fits to the line emission (black curves). The velocity scale is relative to 
the source's redshift of z=5.2979±0.0004, as measured from the molecular line emission. Right: Same, but showing the PdBI data recorded 
with the WideX correlator (CO J=6^5 is re- binned to 20 MHz). 
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=(5.84±0.78), L' 



1) - v „._. , ^ CO(5 _ 4) =(3.70±0.37), 

= (3.82±0.54)xl0 10 K kms^pc 2 . This 

also corresponds to CO J=6— )>5/5— >4 and CO 
J=5— >4/2— >1 line brightness temperature ratios of 
r 65 =1.03±0.16 (i.e., consistent with thermalized) and 
r 52 =0.63±0.10. 

4. ANALYSIS AND DISCUSSION 
4.1. Origin of the CO and Submillimeter Emission 

Within the size of the CO beam, there are three galax- 
ies with colors that are consistent with the CO redshift 
of AzTEC-3. By combining the integrated emission of 
the CO J=5— >4 and 6—^5 lines detected at 12a signifi- 
cance and fitting the u — v data with a circular Gaussian, 
we find that the source is unresolved down to a FWHM 
diameter of 1.0"±0.7" (~6±4kpc at 2=5.3). This is con- 
sistent with a limit of 1.3"j^3„ (~8jl| kpc) as derived 
from the CO J=2— >1 data. Only one of the three color- 
selected galaxies falls within this smaller area, identifying 
the z^26 galaxy targeted by the optical spectroscopy as 
the real counterpart of the CO emission (Fig. [3]). The 
astrometric accuracy of the CO(J=2— >1) and combined 
CO J=5->4 & 6->5 detections is 0.15" and 0.12". The 
peak position of the CO emission is consistent with that 
of the 890 /im contin uum e mission (astrometric accuracy: 
0.14"; Younger et al. [2001 and the 3.6 llth em ission (as- 
trometric accuracy: 0.2"; Sanders et al. I2007h within the 
relative uncertainties (Fig. 13]). Interestingly, the emission 
at all these wavelengths is offset from the HST/ ACS i- 
band image of the source by ^0.5" (3 kpc), but is con- 
sistent with a marginal, much fainter peak in the rest- 



frame UV image (which we interpret to be part of the 
same galaxy; inset in Fig. Given the 0.1" relative 
astrometric accuracy of the HST data, this offset is for- 
mally significant. We conclude that this offset is due 
to a combination of lacking surface brightness sensitiv- 
ity and dust obscuration in the rest-frame UV image, a 
similar effec t as s een in the z=4.055 SMC GN20 (e.g., 
Carilli et al. I2010h . This effect also explains the substan- 
tial difference between the UV- and far- IR-de rived star 
formation rates of AzTEC-3 (Capak et al. l2010l ). suggest- 
ing that the regions of most intense star formation are 
highly dust-obscured. 

4.2. CO Line Excitation 

Based on the CO excitation ladder of AzTEC-3, we can 
constrain the line radiative transfer through Large Veloc- 
ity Gradient (LVG) models, treating the gas kinetic tem- 
perature and density as free parameters. For all calcu- 
lations, the H2 ortho-to-para ratio was fixed to 3:1, the 
cosmic microwave background tempe rature was fixed to 
17.16 K (at z=5.298), and the Flower (|200l CO collision 
rates were used. We adopted a CO abundance per veloc- 
ity gradient of \ CO]/(dv/dr ) = 1 x 10~ 5 p c (km s' 1 )' 1 
(e.g., Weifi et al. [2001 l2007t Riechers et al. [2006b . 

The data are poorly fit by single-component mod- 
els, which underpredict the CO( J=2— sT) flux by ^30%. 
The data can be fitted reasonably well with two gas 
components, which are represented by a 'diffuse', low- 
excitation component with a kinetic temperature of 
7kin=30K and a gas density of p gas =10 2 ' 5 cm -3 , and a 
more 'dense', high-excitation component with Tki n =45K 
and / o gas =10 4,5 cm' 3 in our model (Fig. [4j. Assuming 
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Fig. 3. — Overlays of the CO(J=5— >4) emission (gray co ntours ) 
and SMA 890 fim (left; 2.7" X 2.1" resolution; Younger et al. l2007f ). 
as well as Spi tzer/I RAC 3.6 Atm (right; PSF FWHM = 1.6"; from 
Sanders et al. [2003) emission (black contours in each panel), with 
an HST/ A CS i-band image (gray scale and inset; from Scoville 
etal. l2007l ). 

that the molecular gas is distributed in a face-on, cir- 
cular disk, the observed limits of 1.3" and 1.0" on the 
sizes of these low- and high-excitation gas components 
yield CO disk filling factors of >75% and >10%. The 
'diffuse' gas component has physical properties compara- 
ble to those of nearby spiral galaxies and 'nor mal' h igh-z 
star-forming galaxies (e.g., Dannerbauer et al. 2009), and 
contributes ^35% to the model-predicted CO(J=1^0) 
flux. The 'dense' gas component has properties compara- 
ble to ultra-luminous infrared galaxy (ULIRG) nucl ei and 
high- z FIR-luminous quasars (e.g., Riechers et al. 120061 
2009]), contributing ~65% to the CO(J=1^0) flux. 

The model suggests an excitation- 
corrected CO(J=l— >0) line luminos- 
ity of L' co =6.6xl0 10 K kms^pc 2 , and 
L Fm / L' CO =260 L Q (K kms^pc 2 )" 1 . This lumi- 
nosity ratio is comparable to 'typical' z>2 S MGs and 
quas ar hos t galaxies (e.g., Greve etal. [2005; Riechers 
et al. 12003) . 

4.3. Total Molecular Gas Mass and Dynamical Mass 

Given the gas excitation conditions in AzTEC-3, we 
derive the total molecular gas mass based on a ULIRG 
conversion factor of aco=0.8 Af© (K k ms" 1 p c 2 ) -1 from 
L' co to M gas (Downes & Solomon 1998), yielding 

M gas =5.3xl0 10 M Q O This corresponds to 5.3xM* in 
this system, implying that the baryonic mass is domi- 
nated by the gaseous component. This means that the 
intense, 1800 MQyr -1 starburst in this galaxy has suf- 
ficient fuel to more than sextuple M* throughout the 
remainder of its duration. This also sets the gas deple- 
tion timescale Td e p=Af gas /SFR to 30 M yr, com parable to 
lower reds hift SMGs (e.g., Greve et al. I200"5t Schinnerer 
et al. |2008[) and z>4 quasar host galaxies (e.g., Riechers 
etal.[20(p. 

Adopting the size limit of 1.0" (~6.2kpc), this also 
corresponds to a limiting average gas surface density of 
E gas >1.7x 10 9 Mq kpc~ 2 . This is comparable to the val- 
ues found for z>2 S MGs (typically 2x 10 9 M Q kpc" 2 ; 
Tacconi etal. 120061 [20085, and about half of the 
peak gas density of the z=6.42 quasa r J1148-I-5251 
(3.5xl0 9 M©kpc~ 2 ; Riechers et al. 120091 ) . Assuming an 

10 Assuming a Milky- Way-like a co =3.5M a (K kms -1 pc 2 ) -1 
for the low-excitation component would increase M gas by a factor 
of 2.2. 
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Fig. 4. — CO excitation ladder (points) and LVG model (lines) 
for AzTEC-3. The model (solid line) consists of two gas com- 
ponents (dashed lines): a Milky- Way-like, low-excitation compo- 
nent with a kinetic temperature of Ti;i n =30K and a gas density of 
Pgas=10 2 ' 5 cm - 3 , and a ULIRG-like, high-excitation component 
with T kin =45K and p gas =10 4 5 cm -3 . 

intrinsic S gas in this range (i.e., ^2.7x 10 9 Mq kpc~ 2 ), 
this may indicate that the gas reservoir in AzTEC-3 
has an intrinsic circular gas disk-equivalent radius of 
ro^2.5 kpc. 

Based on the velocity and size information extracted 
from the CO data, we can determine a value for 
the dynamical mass (A/d yn ) of AzTEC-3. Given 
that the contributions from the central supermas- 
sive black hole (M B h<0.1% M* cx 10 7 M q ; Alexan- 
der etal. I2008D and dust (M dust <l-2% M gas ; e.g., 
Michalowski et al. I2010D to the total mass bud- 
get are expected to be small, we here approximate 
M to t~M dyn ~M gas +A4+Af DM , where A/ DM is the con- 
tribution from dark matter. We further define the gas 
mass fraction as fgas—M sas /Md yn and the baryonic mass 
fraction as /bary=l-/DM=[Af gas +M*]/Af dyn . 

In Fig. [51 / gas and /bary are shown as a function of 
ro and the inclination i of the gas disk. For reference, 
r =2.5kpc corresponds to A/d yn sin 2 i=1.4x 10 11 Af , i.e., 
/gas~0.4 and /bary~0.45E3 From the observations of the 
gas and dust in this galaxy, we can infer three main con- 
straints on Mdyn (shaded regions in Fig. [5j. First, the 
CO observations place an upper limit of 1.3" on r$. Sec- 
ond, of the >20 SMGs detected in CO emission to date, 
none has a FWHM linewi dth of >1200 kms" 1 (median: 
530 kms" 1 ; Coppin etal. 2008). Taking this as an up- 
per limit on the intrinsic, inclination-corrected linewidth 
w co=dfFWHM sin -1 ? yields a lower limit of ?>24° (i=0° 
corresponds to face-on). And third, assuming that the 
starburst disk is supported by radiation pressure on 
dust grains yields an Eddington limit on t he FIR flux 
of F pm^l O 13 L Q kpc" 2 (e.g., Scoville HH Thompson 
et al. [2005) . Assuming that all Lfir is due to star for- 
mation in the gas-rich reservoir, this yields a lower limit 

Assuming that the system is seen close to edge-on. 
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Fig. 5. — Constraints on the gas (/gas) and baryonic (/b ary ) 
mass fractions for different line-of-sight inclination angles i for 
AzTEC-3 (where i=90° corresponds to edge-on), ro is the equiv- 
alent radius of a circular disk used for deriving M^ yn . The upper 
limit derived from the CO data and the lower limit derived from 
the Eddington limit on Lpir (fEldd) are indicated by the vertical 
lines and shaded areas. A lower limit on i is determined by as- 
suming that «co=<i''FWHM sin -1 i<1200 kms -1 , as indicated by 
the lower dotted curve and shaded area. The allowed range of i 
for an assumed dark matter fraction of /dm=1 _ /bary=25%±10% 
(corresponding to / gas =63%±8%) is indicated by the thin dotted 
horizontal/arrow-shaped lines. 



on ro of rEdd=735pc. Note that the corresponding dust 
brightness temperature limit of Tjjjj:=88 K is more than 
twice as high as the Td us t obtained fro m fitt ing the spec- 
tral energy distribution (Capak etal. I2010T ) . which may 
indicate that i*FiR in AzTEC-3 is well below the Ed- 
dington limit on average, contrary to what is seen, e.g., 
in th e cent er of the z=6. 42 qua sar J1148+5251 (Walter 
et al. [20091 Riechers et al. [2009h . 

Based on these constraints on Md yn , and assuming 
that aco and M+ are correct, another interesting limit 
arises: for the galaxy not to be dark matter-dominated 
(/bary<50%) within its central few kpc, the system has 
to be seen at i>35°, and the starburst disk has to be rela- 
tively compact (ro<2.3kpc). Assuming a p robab ly more 
realistic / D m=25%±10% (e.g., Daddi et al. [2010h would 
suggest limiting i>44°±4° (at r E dd) and r <1.5±0.2kpc 
(at i=90°). 

Overall, the observations thus may favor a relatively 
compact, highly inclined galaxy, with a high, per- 
haps dominant, fraction of molecular gas (/ gas =0.3-0.8) 
that also dominates the baryonic mass in this system 
(/gas— 0.84 /bary)- Despite the limited spatial resolution 
of this detection experiment, the diagnostic plot intro- 
duced here allows us to constrain the dynamical mass in 
this system to within a factor of *~2-3. However, dynam- 
ically resolved CO observations at high (<0.7") spatial 
resolution are necessary to determine to what degree the 
assumptions made here are correct. Also, potential ef- 
fects of non-virial dynamics (such as in a merger) require 
further investigation. 

5. CONCLUSIONS 

We have detected a molecular gas reservoir in the 
highest-redshift SMG AzTEC-3 (z=5.298) that is com- 
parable in mass t o 'typ ical' z>2 SMGs (e.g., Greve etal. 
120051 Ivison et al. I201C* Harris et al. I2010D and likely has 



a low-excitation component similaiP^l in gas prope rties to 
that seen in another z>4 SMG (Carilli et al. l2010l ). but in 
addition has a high-excitation gas component with prop- 
erties similar_to z>4 quasar host galaxies (e.g., Riechers 
etal. I2006L 120091 ). The peak of the gas distribution is 
coincident with the rest-frame optical and far-infrared 
emission (stellar light and dust-enshrouded star forma- 
tion), but slightly offset from the peak of the rest- frame 
UV emission (unobscured star formation), which appears 
to consist of multiple bright clumps. Given the high gas 
mass fraction and star formation rate of this galaxy, this 
suggests the presence of a heavily obscured starburst, 
possibly triggered by a major, 'wet' (i.e., very gas-rich), 
<8kpc-scale merger. This is reminiscent o f wha t is seen 
in 'typical' z>2 SMGs (e.g., Tacconi et al. 120081 ). 

The detection of luminous CO emission implies rela- 
tively advanced enrichment with heavy elements in the 
material that fuels the observed early burst of star forma- 
tion in AzTEC-3. Assuming a Galactic abundance of CO 
yields A/(CO)~7x 10 7 M Q . Such a level of early metal 
enrichment could be achieved through a few x 10 7 hun- 
dred solar mass population-Ill stars (e.g., Walter etal. 
2003). This ordcr-of-ma gnitudc estimate, however, cor- 
responds to >10% M*, and would require a quite 'top- 
heavy' IMF. We thus speculate that asymptotic giant 
branch stars and supernovae may be extremely efficient 
in enriching their environments at such early epochs, or 
that the stellar mass may be underpredicted due to ob- 
scuration in the most intensely star-forming regions. 

The unusually high-excitation gas component in 
AzTEC-3 (for a SMG) raises the question whether or not 
the environment may play a role in determining its star 
formation properties. AzTEC-3 has two close, massive 
companions with consistent photometric redshifts within 
only 12.2 kpc. The molecular gas reservoir is clearly spa- 
tially separated from the companions; however, gravi- 
tational interactions may still play a role at such small 
distances, which may explain an above average peak gas 
density (leading to high CO excitation). The most sim- 
ilar object probably is the z=4.055 SMG GN20, which 
evolv es in a comparably overdense environment (Daddi 
et al. l2009f ). Interestingly, the CO reservoir in AzTEC-3 
is more compact, (partially) higher excited, and less than 
half as massive as that in GN20. This may either indi- 
cate that the (overdense) cosmic environments of both 
galaxies are considerably different (GN20 even has two 
massive SMG companions within 180 kpc distance), or 
that they are in different stages of their evolution as a 
SMG. The likely presence of a substantial low-excitation 
gas component besides highly-excited gas in both sys- 
tems may point at a fundamental difference in physi- 
cal properties between the star-forming environments in 
z>4 SMGs and z>4 FIR- luminous quasar host galaxies 
(e.g., Riechers et al. I2006h . which may trace overall dif- 
ferent evolutionary stages of galaxies with comparable, 
high gas masses and star formation rates. This differ- 
ence may arise due to a higher average concentration of 
the gas in the quasar hosts (rather than large-scale in- 
fluence of the active nucleus), e.g., due to later merger 
stages. This would be consistent with a scenario where 
SMGs evolve into FIR- luminous quasars (e.g., Sanders 

12 Future CO( J=l— >0) observations are desirable to characterize 
this component in more detail. 
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et al. 119881 : Coppin et al. 12008ft . Both populations are 
likely progenitors of the most massive galaxies in the 
present-day universe. 

This investigation clearly motivates more detailed 
studies of the first SMGs and quasars, which will enable 
us to directly probe the scales that are critical to unravel 
the physical processes that drive the clustered evolution 
of massive galaxies in the early universe, back to the first 
billion years after the Big Bang. 
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